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Abstract 

In  a  numerical/theoretical  effort,  the  dynamics  of  coherent  structures  in  unforced,  and 
periodically  forced,  turbulent  weak  wall  jets  has  been  investigated.  The  objective  of 
this  research  is  to  better  understand  the  effect  of  pulsed  blowing  used  for  separation 
control  on  airfoils.  The  focus  has  been  on  flat  plate  geometries.  In  the  simulations,  a 
weak  wall  jet  is  generated  by  steady  and/or  pulsed  wall-tangential  blowing  through  a 
wall-slot  into  a  constant  pressure  or  adverse  pressure  gradient  boundary  layer.  Numer¬ 
ical  approaches  taken  are  Direct  Numerical  Simulations  (DNS),  Large-Eddy  Simula¬ 
tions  (LES),  and  Unsteady  Reynolds- Averaged  Navier-Stokes  calculations  (URANS). 
Using  three-dimensional  LES,  it  was  shown  that  large-amplitude  pulsed  blowing  can 
prevent  flow  separation.  However,  in  order  to  draw  definite  quantitative  conclusions 
regarding  the  evolution  of  two-dimensional  coherent  structures  generated  by  pulsed 
blowing,  additional  simulations  with  better  grid  resolution  would  be  required.  Such 
simulations  are  computationally  very  expensive.  Therefore,  at  much  lower  computa¬ 
tional  expense,  two-dimensional  URANS  was  employed  for  parameter  studies  of  the 
nonlinear  instability  of  the  weak  wall  jet.  It  was  found  that  forcing  at  high  frequen¬ 
cies  produces  disturbances  that  grow  initially,  while  disturbances  are  always  damped 
for  low  frequency  forcing.  For  large  forcing  amplitudes,  however,  it  appears  that  low 
frequency  disturbances,  while  damped,  persist  over  a  longer  streamwise  distance. 
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1  Motivation  /  Background 

Experimental  research  by  Wygnanski  and  co-workers  (Wygnanski  1997)  has  provided 
considerable  evidence  that  periodically  forced  wall  jets  (tangential  oscillatory  blowing) 
may  be  an  effective  method  for  the  control  of  separation  of  flows  over  single-element 
or  segmented  airfoils  (including  flaps).  However,  many  of  the  fundamental  mecha¬ 
nisms  that  are  responsible  for  the  often  striking  effect  of  oscillatory  blowing  were  not 
understood.  Therefore,  in  close  collaboration  with  the  experimental  effort,  direct  nu¬ 
merical  simulations  (DNS),  large-eddy  simulations  (LES),  and  stability  investigations 
have  being  carried  out  for  typical  wall  jet  configurations.  These  configurations  arise 
when  periodic  forcing  is  applied  for  separation  control  in  practical  applications,  e.g., 
separation  control  for  flows  over  airfoil/flaps.  To  this  end,  we  have  performed  numer¬ 
ical  simulations  for  weak  wall  jets,  that  is,  for  wall  jets  that  are  embedded  in  a  strong 
external  stream.  Future  simulations  are  also  planned  for  wall  jets  over  curved  sur¬ 
faces.  The  work  described  is  a  continuation  of  a  research  effort  funded  previously  by 
AFOSR,  where  the  major  emphasis  has  been  on  strong  wall  jets  (no  external  stream 
or  a  weak  external  stream)  and  for  flat  plate  surfaces  only.  Most  of  the  necessary 
numerical  tools  required  for  the  present  (and  future)  numerical  simulations  were  devel¬ 
oped  during  this  previous  research  effort  (e.g.  Navier-Stokes  codes  for  DNS  and  LES). 

--  -  ) 

The  goal  of  the  present  research  effort  has  been,  with  the  joint  experimental/numerical 
approach,  to  uncover  the  governing  physical  mechanisms  responsible  for  the  striking 
effect  of  oscillatory  blowing  on  separation,  and  thus  provide  the  physical  understand¬ 
ing  required  to  further  develop  this  technique  for  practical  use.  Towards  this  end  we 
have  performed  simulations  using  simplified  model  geometries  that  allow  a  reliable 
investigation  of  the  relevant  mechanisms  and  prevent  contamination  by  effects  that 
are  not  relevant. 
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a)  Experiment  by  Seifert  et  al.  (1994)  b)  Numerical  Simulation  (DNS/LES) 

Figure  1:  Flow  over  a  slotted  flap  with  oscillatory  forcing. 

Over  the  past  decade,  a  great  amount  of  experimental  work  has  been  performed  by 
Wvgnanski  and  coworkers  (Wygnanski  1997)  which  clearly  demonstrates  the  effec¬ 
tiveness  of  using  pulsed  wall  jets  for  delaying  separation.  The  setup  of  one  such 
experiment  (Seifert  et  al.  1994)  is  illustrated  in  Figure  la.  In  the  experiment,  a  high 
momentum  jet  that  is  blown  through  the  slot  along  the  upper  surface  of  the  flap 
is  periodically  forced  at  the  location  indicated  in  Figure  la.  The  additional  forcing 
strongly  improves  the  effectiveness  of  the  weak  wall  jet  in  preventing  separation,  and 
thus  allows  for  significantly  higher  angles  of  attack.  The  effectiveness  of  using  pulsed 
wall  jets  for  delaying  separation  has  been  clearly  demonstrated  in  experiments  (Seifert 
et  al.1993,  1996,  Wygnanski  1997). 

While  use  of  LES  and  DNS  for  flows  over  actual  slotted  flaps  are  the  final  goal  of  our 
computational  research  (Figure  lb),  several  intermediate  steps  must  be  taken  first  to 
explore  different  aspects  of  this  complex  flow  geometry.  Clearly,  for  wall  jets  over 
actual  slotted  flaps,  the  combined  effects  of  adverse  pressure  gradient  and  curvature 
plays  a  major  role.  To  isolate  the  relevant  mechanisms,  this  complex  flow  geometry 
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Figure  2:  Generic  wall  jet  configurations  for  investigating  different  aspects  of  the  flow 
over  a  slotted  flap. 


is  first  broken  down  into  simpler  modules  so  that  the  effect  of  pressure  gradient  and 
curvature  can  be  investigated  separately. 


These  intermediate  steps  are  outlined  in  Figure  2.  During  the  first  year  of  the  present 
research  project  we  have  focussed  on  turbulent  weak  wall  jets  with  zero  free  stream 
pressure  gradient  (Figure  2a).  Typical  results  from  these  simulations  will  be  discussed 
below.  In  the  second  year,  we  have  investigated  turbulent  weak  wall  jets  subjected  to 
strong  adverse  pressure  gradients  (Figure  2b).  The  effect  of  curvature  will  be  studied 
as  part  of  the  new  research  project  (AFOSR  GRANT  F49620-00-1-0069),  at  first  for  a 
strong  wall  jet  over  a  curved  wall  geometry  (constant  curvature,  Coanda  flow  (Figure 
2c).  Ultimately,  the  entire  flow  over  a  slotted  flap  will  be  investigated  where  the 
effects  of  pressure  gradient  and  curvature  are  present  simultaneously  (Figure  2d). 
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2  Computational  Method 

The  Navier-Stokes  codes  employed  for  our  numerical  investigations  were  originally 
developed  in  our  research  group  for  boundary  layer  transition  (see  Meitz  and  Fasel, 
2000)  and  have  since  been  adapted  to  the  wall  jet  flow  geometry. 

2.1  Governing  Equations 

In  our  codes,  the  complete  incompressible  Navier-Stokes  equations  are  solved  in 
vorticity-velocity  formulation.  Regardless  of  which  form  of  time  averaging  or  spa¬ 
tial  filtering  is  applied,  the  resulting  equations  take  on  the  form, 

f  +  vx(5xC-v.?)  =  ^.  (i) 

The  overbar  on  the  velocity  vector  V  and  on  the  vorticity  vector  u  denotes  the 
operation  of  averaging  the  appropriate  components  of  the  corresponding  vector.  The 
quantity  t  =  [r^j  is  the  turbulent  stress  tensor  which  must  be  modeled.  Equation  (1) 
is  non-dimensionalized  by  the  free-stream  velocity  Uoo,  a  reference  length  L,  and  a 
global  Reynolds  number  defined  as  Re  =  UaoL/i>  where  v  is  the  kinematic  viscosity. 
Tu  .addition  to  the  transport  equations  for  the  vorticities,  the  governing  equations 
_  also  include  a  set  of  Poisson  equations  for  the  velocities.  From  the  definition  of 
vorticity  and  the  fact  that  both  the  velocity  and  the  vorticity  fields  are  solenoidal, 
three  equations  for  the  velocity  components  are  obtained  (for  more  details  see  Meitz 
and  Fasel,  2000).  The  transport  equations  for  the  vorticities  and  the  Poisson  equations 
for  the  velocities  form  a  nearly  complete  system  of  governing  equations.  In  order  to 
obtain  closure,  a  set  of  relations  is  needed  for  the  turbulent  stress  tensor  [r^].  These 
equations  are  provided  in  the  form  of  our  new  Flow  Simulation  Methodology  (FSM), 
see  section  2.3. 


6 


2.2  Smagorinsky  Subgrid-scale  Model 

For  the  LES  of  turbulent  wall  jets,  as  a  first  step,  a  basic  Smagorinsky-type  subgrid- 
scale  model  (suggested  by  Speziale  1995)  was  implemented  into  the  Navier-Stokes  code 
using  a  constant  Smagorinsky  coefficient.  The  model  assumes  that  the  sgs  stresses 
follow  a  gradient-diffusion  process,  similar  to  molecular  motion.  Thus, 


Tij  =  T^ES  =  -2 vrSij, 

(2) 

where  vT  is  the  effective  viscosity  of  the  subgrid  scales, 

»t  =  P'JzSijSij, 

(3) 

where 

l  =  Cs(l-  e~(y+/A+^  {AxAyAz)1/3 

(4) 

and  Cs  =  0.065  is  the  Smagorinsky  constant.  Ax,  Ay  and  Az  are 

the  grid  widths  in 

the  x ,  y  and  2  directions,  respectively. 

This  baseline  eddy-viscosity  model  performs  reasonably  well  for  LES  of  zero-pressure 
gradient  weak  wall  jets.  However,  from  preliminary  calculations  we  know  that  this 
baseline  model  cannot  properly  estimate  the  subgrid-scale  stress  for  more  complex 
flows  (such  as  wall  jet  flows  with  strong  adverse  pressure  gradients)  which  may  be 
close  to  separation  (Figure  2c).  Therefore,  we  have  implemented  a  new  subgrid-scale 
model  into  the  wall  jet  code. 
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2.3  New  Flow  Simulation  Methodology  (FSM) 

In  collaboration  with  C.  Speziale  (1998),  we  have  been  developing  a  new  Flow  Simula¬ 
tion  Methodology  (FSM)  to  be  applicable  for  complex  turbulent  flows.  A  key  feature 
of  this  new  FSM  is  that  it  is  consistent  with  DNS  for  fine  grid  resolution  and  with 
Reynolds-Averaged  Navier-Stokes  (RANS)  calculations  for  coarse  grid  resolution. 


This  is  achieved  by  modeling  of  the  turbulent  stress  tensor  [r^]  used  in  equation  (1) 
as 

Tii  =  m /ljt?,  (5) 

where  t-j  is  the  Reynolds  stress  tensor,  and  f(A/Lk)  is  the  so-called  contribution 
function.  The  terms  A  =  [(dx2  +  dy 2  +  dz2)/3]1/2  and  Lk  =  ^4/el/4  are  the  repre¬ 
sentative  computational  grid  size  and  the  Kolmogorov  length  scale,  respectively.  In 
the  coarse  grid  limit,  f(A/Lk)  -+  1,  a  RANS  calculation  is  recovered.  In  the  fine  grid 
limit,  f((A/Lk)  — >  0,  the  FSM  approaches  a  Direct  Numerical  Simulation  (DNS).  A 
detailed  discussion  of  how  to  design  an  appropriate  contribution  function  f(A/Lk)  is 
given  by  Zhang  et  al.  (2000). 


For  the  Reynolds  stress  t*  in  equation  (o),  the  Algebraic  Stress  Model  (ASM)  of 
Gatski  and  Speziale  (1993)  is  used  that  allows  for  a  misalignment  between  the  stress 
and  strain  rate  tensor  as  it  occurs  in  non-equilibrium  flows. 


t*  =  2-K6tj  -  /( 7J, 0  a +  t*2~2  {SikWkj  +  SjkWki ) 

(SikSkj  —  -SkiSki5ij 


(6) 


Regularization  is  performed  on  the  function  f(r),£)  to  remove  the  possibility  of  divi¬ 
sion  by  zero  yielding, 

3  _  3(1  +  7 f) 

3  -  2r?2  +  6£2  ~  3  +  T]2  +  6r 2  ' 


/(* 7.0  = 


(7) 
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The  terms  77  and  £  are  the  irrotational  strain  rate  invariant  and  the  rotational  strain 
rate  invariant,  respectively.  They  depend  on  Sij  and  Wij  in  the  following  way: 

(«> 

{  =  —  (WijWij)'12  - .  (9) 

The  values  of  coefficients  appearing  in  equations  (6),  (8)  and  (9)  are  ai  =  0.227, 
a2  =  0.0423,  and  a3  =  0.0396. 

The  quantities  k  and  e  appearing  in  equation  (6)  are  the  turbulent  kinetic  energy  and 
the  turbulent  dissipation  rate.  They  are  computed  from  a  standard  set  of  turbulent 
transport  equations, 


Auxiliary  relations,  including  an  expression  for  the  eddy-viscosity  i *r  and  wall  damp¬ 
ing  functions  are, 

=  /e,  U  =  l-e-»+'A+,  (12) 

/e2  =  l  _  e-s/Wfc/io  (13) 

The  constants  needed  for  the  computation  of  equations  (10)  through  (12)  are, 

C„  =  0.09,  C(l  =  1.44,  Ct2  =  1.83, 
ok  =  1,  o(  =  1.3,  A+  =  25. 
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2.4  Numerical  Scheme 

A  detailed  review  of  the  numerical  schemes  used  for  the  solution  of  equation  (1) 
and  the  Poisson  equations  for  the  velocities  can  be  found  in  Meitz  and  Fasel  (2000). 
In  summary,  the  time  dependent  spatially  varying  flow  field  is  calculated  using  the 
so-called  spatial  model  (inflow-outflow).  For  the  time  integration  a  fourth-order  accu¬ 
rate  explicit  Runge-Kutta  scheme  is  employed.  Spatial  derivatives  in  the  streamwise 
direction  are  approximated  by  fourth-order  compact  differences  employing  upwind 
biased  and  downwind  biased  switching.  Wall-normal  derivatives  are  computed  with 
fourth-order  compact  differences  that  allow  for  grid  stretching.  Pseudo-spectral  de¬ 
composition  is  utilized  in  the  spanwise  direction,  z.  Lastly,  a  fast  Poisson  solver  is 
implemented  for  solving  the  Poisson  equations  for  the  velocity  components. 

The  FSM  requires  simultaneous  solution  of  the  k-e  equations  (see  section  2.3).  Space 
and  time  dependent  values  for  k,e  are  required  for  the  contribution  function  /(A/ Lk) 
(to  determine  Lk )  and  for  modeling  of  r-j.  The  time  dependent  k  and  e  equations 
(see  section  2.3)  are  solved  using  an  ADI  method  with  second-order  accuracy  both 
in  space  and  in  time.  The  spatial  derivatives  of  the  turbulent  stresses  are  computed 
using  the  same  difference  approximation  as  for  the  computation  of  the  resolved  scales 
as  discussed  above. 
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3  Validation  of  the  new  FSM 

The  Navies-Stokes  code  used  for  all  our  computations  has  been  extensively  validated 
in  the  past  (Meitz  and  Fasel,  2000).  Also,  the  code  used  for  our  Large-Eddy  Simula¬ 
tions  using  the  Smagorinsky  model  has  been  tested  for  numerous  flow  configurations 
(  see  e.g.  Fasel  et  al.,  1997).  However,  the  code  extensions  required  for  our  URANS 
calculations  and  the  new  FSM  needed  to  undergo  a  series  of  benchmark  calculations. 
In  these  calculations,  the  code’s  capability  to  capture  correctly  the  various  aspects  of 
the  present  very  complex  flow  configuration  was  thoroughly  tested. 

3.1  Turbulent  Boundary  Layer 

The  performance  of  the  FSM  code  has  been  evaluated  extensively  for  turbulent  bound¬ 
ary  layers  without  and  with  adverse  pressure  gradient  in  the  RANS  limit  and  using 
the  FSM  as  an  LES.  Results  from  these  investigations  have  been  published  by  Zhang 
et  al.  (2000a,  2000b).  Both,  mean  flow  and  turbulence  statistics,  match  experimental 
findings  reported  in  the  literature  very  well.  One  test  case  for  FSM  calculations  in  the 
RANS  limit  is  now  briefly  discussed,  since  the  turbulent  boundary  layer  computed  in 
these  particular  calculations  also  serves  as  the  base  flow  for  the  URANS  calculations 
m  section  5.  In  Figure  3,  profiles  for  several  flow  quantities  from  RANS  calculations 
with  various  turbulence  models  are  compared  in  wall-coordinates  {Reg  =  2525).  The 
turbulence  models  examined  are  the  Baldwin-Lomax  model,  two  standard  k  —  t  mod¬ 
els  using  near-wall  corrections,  MK  and  SSA  [Myong  and  Kasagi  (1990)  and  Sarkar 
and  So  (1997)],  and  the  ASM  model  by  Gatski  and  Speziale  (1993). 1  As  shown  in  Fig¬ 
ure  3a,  for  all  turbulence  models  considered  the  streamwise  velocities  agree  well  with 
the  log-law.  All  models  are  also  in  close  agreement  for  the  Reynolds  stress  (Figure 
3b),  except  very  close  to  the  wall  {y+  <  20).  The  discrepancy  near  the  wall  is  much 

1  Calculations  using  the  boundary  layer  approximation  (ASM-RABL)  and  using  the  full  Navier- 
Stokes  equations  (ASM-RANS). 
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more  pronounced  for  the  turbulent  kinetic  energy  k  and  the  turbulent  dissipation 
rate  t  shown  in  Figures  3b  and  3c  respectively.  For  comparison,  results  from  a  DNS 
by  Spalart  (1988)  and  from  experimental  measurements  by  Karlsson  and  Johansson 
(1988)  are  also  shown  in  these  graphs.  For  y+  <  100,  the  models  by  Myong  and  Kasagi 
and  Sarkar  and  So  clearly  agree  much  better  with  the  DNS  and  the  experiment.  This 
has  to  be  expected,  since  both  models  employ  empirical  near-wall  damping  functions 
that  are  tuned  to  match  the  asymptotic  near- wall  behavior  of  the  DNS.  The  ASM, 
on  the  other  hand,  minimizes  the  use  of  wall-damping  functions  at  the  price  of  an 
incorrect  asymptotic  behavior  near  the  wall.  However,  in  all  our  calculations  this  de¬ 
ficiency  had  no  impact  on  the  quality  of  the  flow  field,  clearly  because  the  Reynolds 
stress  computed  using  k  and  e  is  very  small  where  the  turbulence  models  disagree 
significantly  (see  close  up  in  Figure  3b).  In  summary,  the  ASM  used  in  our  URANS 
calculation  performs  very  well  for  the  turbulent  boundary  layer  without  the  use  of 
empirical  near-wall  corrections. 
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a)  streamwise  velocity 


c)  turbulent  kinetic  energy 


b)  Reynolds  stress 


d)  turbulent  dissipation  rate 


Figure  3:  Performance  of  the  RANS-solver  for  the  turbulent  boundary  layer  without 
wall  blowing.  Results  from  several  RANS  calculations  using  different  turbulence  mod¬ 
els  are  compared.  Also  shown  are  turbulent  kinetic  energy  and  turbulent  dissipation 
rate  from  a  DNS  by  Spalart  (1988)  and  from  experimental  measurements  by  Karlsson 
and  Johansson  (1988). 
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3.2  Turbulent  Strong  Wall  Jet 

The  turbulent  strong  wall  jet  represents  another  crucial  test  case,  since,  in  addition  to 
a  boundary  layer,  it  includes  one  free  shear  layer  (as  opposed  to  two  free  shear  layers 
for  the  weak  wall  jet).  Unlike  for  weak  wall  jets,  there  are  numerous  experimental  and 
computational  investigations  available  for  comparison.  As  an  example,  a  comparison 
with  the  experiments  by  Eriksson  et  al.  (1998)  for  a  jet  exit  Rej  =  10000  is  shown  in 
Figure  4  for  the  same  k  —  e  models  (MK,  SSA,  and  ASM)  and  the  same  flow  quantities 
as  for  the  turbulent  boundary  layer  discussed  in  the  previous  paragraph.  Again,  ASM 
agrees  well  with  MK  and  SSA  for  streamwise  velocity  (Figure  4a)  and  Reynolds  stress 
(Figure  4b),  despite  the  asymptotically  incorrect  behavior  for  k  and  e  (Figures  4c  and 
d).  In  fact,  all  turbulence  models  exhibit  similar  characteristics  as  for  the  turbulent 
boundary  layer,  except  they  all  undershoot  the  streamwise  velocity  profile  from  the 
experiment  which  is  a  well-known  deficiency  of  two-equation  models. 

In  order  to  demonstrate  the  capability  of  our  method  to  accurately  capture  large 
coherent  structures  in  wall  jet  flows,  we  have  performed  a  detailed  comparison  with 
experiments  by  Katz  et  al.  (1992)  on  forced  strong  wall  jets.  The  results  were  ex¬ 
tremely  encouraging.  Not  only  were  we  able  to  match  mean  flow  and  eigenfunctions 
of  the  experiments  extraordinarily  well,  but  we  could  also  provide  a  possible  explana¬ 
tion  for  the  mismatch  between  the  dominating  measured  frequency  in  the  inner  and 
the  outer  region  of  the  wall  jet.  Our  findings  are  detailed  in  a  paper  presented  at  the 
Fluids  2000  conference  in  Denver,  CO  (Seidel  and  Fasel,  2000).  A  copy  of  the  paper 
is  attached  to  this  report. 
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Figure  4:  Performance  of  the  RANS-solver  for  the  turbulent  strong  wall  jet.  Results 
from  several  RANS  calculations  using  different  turbulence  models  are  compared.  Also 
shown  are  experimental  measurements  by  Eriksson  et  al.  (1998). 
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3.3  Turbulent  Weak  Wall  Jet 

For  further  testing  of  the  model  we  performed  RANS  calculations  of  a  turbulent  weak 
wall  jet  which  was  generated  by  blowing  into  a  turbulent  boundary  layer  through  a 
slot  at  a  shallow  angle  (geometry  in  Figure  2b).  Over  the  blowing  slot,  streamwise 
and  wall-normal  velocity  profiles  are  specified.  This  simple  approach  for  modeling 
the  blowing  slot,  while  allowing  for  some  flow  adjustment  close  to  the  slot,  produces 
the  correct  downstream  development  of  the  wall  jet,  as  was  demonstrated  in  earlier 
simulations  (DNS  of  laminar  wall  jets,  LES  of  turbulent  wall  jets  using  the  Smagorin- 
sky  baseline  model). 

For  the  new  LES,  a  Dirichlet  boundary  condition  is  also  imposed  for  the  turbulent 
kinetic  energy  ( ksiot  >  0).  ksiot  provides  a  measure  for  the  turbulence  level  within  the 
injected  fluid.  We  have  verified  that  the  choice  of  ksiot  only  influences  the  flow  close 
to  the  blowing  slot  and  is  not  crucial  for  the  development  of  the  weak  wall  jet  further 
downstream.  This  is  illustrated  in  Figure  5,  where  the  streamwise  development  of 
the  wall  jet  velocity  is  compared  for  three  different  RANS  calculations:  ASM  with 
ksiot  =  0.012  (case  1),  standard  k  —  e  model  with  ksiot  =  0.012  (case  2),  standard  k  —  e 
model  with  ksi0t  =  0.002  (case  3).  While  the  three  cases  differ  close  to  the  blowing 
slot,  they  are  all  in  good  agreement  further  downstream.  Over  the  blowing  slot,  the 
standard  k  —  e  model  is  certainly  not  well  calibrated,  and,  as  a  result,  it  is  very  dis¬ 
sipative  (case  2).  Consequently,  the  weak  wall  jet  initially  decays  more  rapidly  than 
for  the  ASM  (case  1).  This  can  be  offset  by  reducing  ksiot.  The  k  —  e  model  (case  3) 
then  matches  the  ASM  result  (case  1)  even  close  to  the  slot. 

An  important  feature  of  the  new  FSM  is  the  fact  that  it  allows  for  unsteady  calcu¬ 
lations  even  in  the  RANS  limit.  This  is  accomplished  by  the  use  of  a  time-accurate 
Navier-Stokes  code  and  the  use  of  a  self-adjusting  turbulence  model  (ASM).  An  ex¬ 
ample  is  provided  in  Figure  6  which  shows  perspective  plots  of  the  streamwise  velocity 
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Figure  5:  Test  case  for  the  new  LES  methodology:  RANS  calculation  of  a  weak  wall 
jet  blown  into  a  turbulent  boundary  layer  through  a  slot  at  a  shallow  angle  (5.4°). 
Comparison  of  the  streamwise  velocity  components  computed  with  ASM  and  with  a 
standard  k  —  e  model. 


a)  ASM,  kslot=0.012  b)  ASM,  ksiopO.005 


Figure  6:  Unsteady  RANS  calculation  of  weak  wall  jet  using  ASM.  Shown  are  per¬ 
spective  plots  of  the  streamwise  velocity  for  two  different  levels  of  turbulent  kinetic 
energy  that  are  specified  over  the  blowing  slot  ( ksi0t )■ 
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for  two  RANS  calculation  using  ASM  with  different  levels  of  ksiot.  While  in  the  case 
of  kaiot  =  0.012  (Figure  6a)  the  flow  converges  to  a  steady  state,  in  the  the  case  of 
ksiot  =  0.002  (Figure  6b),  downstream  traveling,  unsteady  structures  are  observed. 
These  structures  originate  in  a  small  separation  bubble  which  has  formed  right  down¬ 
stream  of  the  blowing  slot.  Note,  that  for  smaller  ksiot  the  jet  has  less  turbulent 
intensity  and,  as  a  consequence,  turbulent  mixing  is  decreased.  Therefore,  locally 
the  jet  can  separate  from  the  surface.  When  the  standard  k  —  e  is  used  instead  of 
the  ASM,  a  separation  bubble  develops  only  for  a  much  smaller  ksiot,  and  even  when 
vortical  structures  are  eventually  generated,  they  decay  very  rapidly. 


4  LES  using  the  Smagorinsky  Model 


In  parallel  with  the  development  of  our  code  employing  the  new  LES  methodology,  we 
have  investigated  turbulent  weak  wall  jets  using  LES  with  the  Smagorinsky  baseline 
model.  While  these  simulations  certainly  could  not  capture  all  aspects  of  the  complex 
flow  correctly  in  a  quantitative  sense,  they  provided  us  with  valuable  reference  data 
for  our  FSM  approach. 


4.1  Zero  Pressure  Gradient  Flow 

The  weak  wall  jets  were  generated  by  blowing  through  an  angled  slot  in  the  wall  into 
a  developed  flat-plate  boundary  layer  with  zero  free  stream  pressure  gradient.  The 
boundary  layer  in  these  simulations  closely  matches  the  artificially  thickened  bound¬ 
ary  layer  in  the  experimental  investigations  of  forced  turbulent  weak  wall  jets  by 
Wygnanski  and  coworkers  (see  Weidemann,  1996).  An  impression  of  the  overall  flow- 
field  computed  with  our  simulations  is  obtained  from  Figure  7,  in  which  greyscales  of 
instantaneous  spanwise  vorticity  are  displayed.  In  the  simulation,  a  laminar  boundary 
layer  is  tripped  by  3D  random  blowing  and  suction  through  a  slot  close  to  the  inflow 
(Figure  7).  Further  downstream,  the  weak  wall  jet  is  generated  by  tangential  blowing 


1Z9  I A 
3D  random  j  I 
blowing  &  suction 


2D  wall-tangential 
blowing  (5.4°) 


13,5  13.6 

x-direction,  [m] 


Figure  7:  LES  of  a  weak  wall  jet  generated  by  tangential  blowing  through  a  slot  in 
the  wall  into  a  turbulent  boundary  layer.  Shown  are  greyscales  of  the  2D  Fourier 
component  of  the  instantaneous  spanwise  vorticity. 
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(angle  5.4°)  into  the  transitioning  boundary  layer  through  a  second  slot.  In  agree¬ 
ment  with  the  experiments  the  jet  persists  for  some  streamwise  distance  25 cm) 

and  then  disintegrates. 

The  streamwise  development  of  this  weak  wall  jet  can  be  nicely  observed  from  the 
time-averaged  streamwise  velocity  component,  as  illustrated  in  Figure  8  using  per¬ 
spective  plotting.  The  velocity  maximum  of  the  wall  jet  decays  rapidly  and  disappears 
about  25cm  downstream  of  the  slot.  For  comparison,  the  boundary  layer  without  a 
wall  jet  is  shown  in  Figure  8a.  In  the  fully  turbulent  regime  close  to  the  outflow, 
the  velocity  profiles  in  Figures  8a  and  8b  are  virtually  identical.  Results  from  other 
simulations,  for  which  the  wall  jet  was  injected  either  farther  upstream  or  farther 
downstream  compared  to  the  case  of  Figure  8b,  are  shown  in  Figures  8c  and  8d. 
These  results  indicate  that  weak  wall  jets  that  are  either  generated  in  the  regime 
of  early  transition  (Figure  8c)  or  in  the  fully  turbulent  regime  (Figure  8d)  persist 
noticeably  longer  («  35cm).  The  faster  decay  for  the  case  in  Figure  8b  is  probably 
due  to  the  fact  that  the  fluctuations  in  the  flow  are  strongest  during  the  later  stages 
of  transition.  The  destructive  effect  of  large  fluctuations  on  a  weak  wall  jet  were  also 
observed  in  earlier  simulations  (Wernz  and  Fasel,  1997)  where  we  attempted  to  gen¬ 
erate  a  turbulent  weak  wall  jet  by  tripping  a  laminar  weak  wall  jet.  As  an  example 
from  these  simulations,  in  Figure  9  corresponding  perspective  plots  of  the  streamwise 
velocity  are  shown  for  the  initial  laminar  weak  wall  jet  and  the  transitional  weak  wall 
jet.  While  the  unforced  laminar  weak  wall  jet  persists  throughout  the  computational 
domain  (Figure  9a),  when  forced  by  3D  random  blowing  and  suction,  the  wall  jet 
decays  rapidly  as  it  transitions  (Figure  9b). 
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c)  weak  wall  jet  generated  farther  upstream  d)  weak  wall  jet  generated  farther  downstream 

Figure  8:  Comparison  of  weak  wall  jets  generated  by  wall-tangential  blowing  through 
a  slot  at  various  downstream  locations.  Shown  are  perspective  plots  of  the  time- 
averaged  streamwise  velocity  component  (steady  mean  flow). 


a)  laminar  weak  wall  jet  b)  transitional  weak  wall  jet 

Figure  9:  Comparison  of  a  laminar  and  a  transitional  weak  wall  jet.  Shown  are 
perspective  plots  of  the  time-averaged  streamwise  velocity  component  (steady  mean 
flow). 
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4.2  Adverse  Pressure  Gradient  Flow 

Towards  the  goal  of  investigating  tangential  blowing  as  a  tool  for  separation  preven¬ 
tion  on  an  airfoil  geometry,  we  have  studied  the  effect  of  tangential  blowing  into  a 
turbulent  boundary  layer  with  an  adverse  pressure  gradient.  A  typical  case  of  such 
investigations  is  presented  in  Figure  10.  In  the  absence  of  blowing,  the  turbulent 
boundary  layer  separates  due  to  the  strong  adverse  pressure  gradient,  as  illustrated 
in  Figure  10a.  When  steady  blowing  is  employed  (Figure  10b),  the  flow  remains  at¬ 
tached,  although  it  gets  close  to  separation.  With  additional  periodic  forcing  (Figure 
10c)  the  flow  does  not  come  close  to  separation.  This  can  be  seen  most  clearly  from 
the  skin-friction  coefficient  plotted  in  Figure  10  for  the  three  cases.  While  forcing 
the  weak  wall  jet  clearly  produces  two-dimensional  coherent  structures  (Figure  10c), 
their  growth  rate  is  lower  than  could  be  expected  from  experimental  findings  (Zhou, 
personal  communication).  We  suspect  that  the  baseline  Smagorinsky  model,  which 
was  employed  in  this  LES,  provided  too  much  dissipation  and  artificially  dampens 
the  coherent  structures.  Also,  with  the  computational  resources  at  our  disposal,  we 
were  unable  to  provide  sufficient  resolution  for  the  near-wall  region. 

Thus,  for  detailed  investigation  of  the  forced  flow,  there  is  a  clear  need  for  an  improved 
modeling  approach  that  does  not  artificially  dampen  out  the  coherent  structures 
in  the  flow  and  is  efficient  enough  to  allow  for  parameter  studies  with  a  moderate 
computational  expense. 
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a)  turbulent  boundary  layer  without  wall  jet 


b)  turbulent  boundary  layer  with  unforced  wall  jet 


+  periodic  forcing  (f=110Hz) 


c)  turbulent  boundary  layer  with  forced  wall  jet 

Figure  10:  Turbulent  flat-plate  boundary  layer  subjected  to  strong  adverse  pressure 
gradient.  Separation  prevention  using  a  weak  wall  jet.  Shown  are  color-contours  for 
the  2-D  Fourier  component  of  the  instantaneous  spanwise  vorticity. 
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Figure  11:  Color-contours  for  the  2-D  Fourier  component  of  the  instantaneous  stream- 
wise  velocity  for  the  three  cases  shown  in  Figure  10. 
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Figure  12:  Streamwise  development  of  time-averaged  flow  quantities  for  the  three 
cases  shown  in  Figure  10. 
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5  RANS  Calculations 

The  objective  of  our  two-dimensional  RANS  calculations  is  to  investigate  the  gener¬ 
ation  and  evolution  of  two-dimensional  coherent  structures  in  the  weak  wall  jet  at  a 
fraction  of  the  computational  and  logistical  cost  of  three-dimensional  LES  and  DNS. 
Results  for  a  zero  pressure  gradient  weak  wall  jet  which  approximately  matches  an 
experiment  by  Weidemann  (1996)  will  now  be  discussed.  The  strategy  for  computing 
turbulent  weak  wall  jets  using  RANS  follows  the  approach  taken  for  our  LES  (see  sec¬ 
tion  4).  A  turbulent  boundary  layer  is  computed  first,  then  the  wall  jet  is  generated 
by  wall-tangential  blowing  through  a  slot.  For  the  RANS  calculations  presented  here, 
the  Reynolds  number  is  Re0  =  2430  at  the  inflow.  2  The  Reynolds  number  at  the 
blowing  slot,  Reg  =  3420,  and  the  wall  blowing  parameters  (steady  blowing  velocity 
As  =  1.3  •  Ufs  =  19.5m/s,  blowing  angle  a  =  5.4°)  match  those  in  the  LES. 

5.1  Steady  RANS  Calculations  of  the  Unforced  Flow 

With  our  RANS  calculations  we  can  meet  the  grid  resolution  requirements  for  ac¬ 
curately  computing  weak  wall  jets.  Results  from  a  resolution  study  for  the  steady 
flow  are  presented  in  Figure  13.  Three  grid  resolutions  are  considered,  the  lowest  for 
base  flow  1,  the  highest  for  base  flow  3.  The  streamwise  stepsize  for  base  flows  2  and 
3  matches  that  of  the  LES,  for  base  flow  1  it  is  five  times  larger  than  for  the  LES. 
In  wall-normal  direction,  base  flow  1  has  the  resolution  used  in  the  LES,  base  flow 
2  has  twice  and  base  flow  3  three  times  that  resolution.  In  Figure  13a,  the  stream- 
wise  velocity  profiles  for  the  three  base  flows  and  for  the  turbulent  boundary  layer 
are  plotted  for  three  downstream  locations.  Overlayed  are  the  locations  of  the  jet 
velocity  maximum  Umax  and  the  velocity  minimum  Um in.  The  jet  extends  to  about 
150mm  downstream  from  the  center  of  the  blowing  slot.  Farther  downstream,  the 

2 In  section  3.1  mean  flow  data  close  to  the  inflow  are  compared  with  experimental  and  DNS 


results. 
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Figure  13:  RANS  calculation  of  steady  weak  wall  jet.  Shown  are  for  three  different 
grid  resolutions:  a)  streamwise  velocity  profiles  at  three  streamwise  locations  plus  the 
wall-normal  location  of  velocity  maximum  and  minimum,  b)  wall-normal  maximum 
and  minimum  of  streamwise  velocity,  c)  skin  friction  coefficient. 
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velocity  maximum  and  minimum  have  disappeared,  but  a  distortion  of  the  boundary 
layer  still  persists.  Figure  13b  displays  the  streamwise  development  of  the  velocity 
maximum  and  minimum  and  Figure  13c  shows  the  skin  friction  coefficient.  For  all 
quantities  shown  in  Figure  13a-c,  base  flows  2  and  3  match  closely,  while  base  flow  1 
deviates  slightly,  particularly  for  the  skin  friction  close  to  the  blowing  slot.  This  is 
because  the  blowing  slot  is  not  well  resolved  with  only  five  grid  points  (see  Figure  15b). 

In  conclusion,  for  the  steady  flow  the  resolution  used  for  base  flow  2  is  the  best  choice, 
the  higher  wall-normal  resolution  for  base  flow  3  is  unnecessary.  While  the  resolution 
for  base  flow  1  is  adequate  for  steady  flow  calculations,  this  is  not  the  case  when 
computing  the  forced  flow  using  URANS  (see  below). 

5.2  Unsteady  RANS  Calculations  of  the  Forced  Flow 

As  a  first  step  for  gaining  insight  into  the  instability  mechanisms  of  a  weak  wall  jet 
we  have  recomputed  a  pulsed  wall  jet  flow  with  URANS  using  the  same  forcing  pa¬ 
rameters  as  for  the  LES  in  section  4.  For  this  case  (denoted  as  case  1L),  the  forcing 
amplitude  is  very  large,  Ap  =  0.65  •  U/s,  which  is  50%  of  the  steady  blowing  arm 
plitude.  A  forcirg  frequency  of  /  =  IIO.8//2  is  used.  The  URANS  calculation  was 
performed  for  the  three  base  flows  discussed  above  in  section  5.1.  Figure  14  illustrates 
the  resulting  flowfield  for  base  flow  1  with  color  contours  of  instantaneous  spanwise 


x-xslot  tmm] 

Figure  14:  URANS  calculation  of  a  pulsed  weak  wall  jet  (base  flow  1,  medium  grid 
resolution).  Shown  are  color  contours  of  instantaneous  spanwise  vorticity  for  case  1L. 
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Figure  15:  URANS  calculation  of  pulsed  weak  wall  jet,  (case  1L —  steady  blowing 
with  As  =  1.3 Ufs,  pulsed  blowing  with  Ap  =  0.65U/S,f  =  110.87/2).  Shown  are  for 
three  grid  resolutions  the  streamwise  mean  flow  profiles  of  the  unforced  and  forced 
flow,  and  the  Fourier  amplitudes  of  the  streamwise  disturbances.  The  profiles  are 
plotted  for  three  streamwise  locations  (a-c)  corresponding  to  those  in  Figure  13a. 


vorticity.  At  this  particular  time  instant;  one  localized  dipolar  structure  can  be  ob¬ 
served  150mm  downstream  of  the  slot.  The  subsequent  structure  at  40mm  is  in  the 
process  of  formation.  Due  to  the  large  forcing  amplitude,  the  resulting  disturbances 
are  highly  nonlinear  with  localized  vortical  structures.  For  base  flows  2  and  3,  the 
streamwise  resolution  provides  ten  grid  points  for  resolving  these  localized  structures, 
while  base  flow  1  only  provides  two,  which  is  clearly  inadequate. 

Consequently,  as  shown  in  Figure  15,  the  fundamental  disturbance  amplitudes  from 
the  URANS  calculation  with  base  flow  1  strongly  deviate  from  those  from  calcula¬ 
tions  with  base  flows  2  and  3.  Figure  15  also  indicates  that  the  disturbances  are 
strongly  damped.  The  disturbance  amplitudes  in  Figure  15c  are  one  order  of  magni¬ 
tude  lower  than  in  Figure  15a.  This  is  not  really  surprising  for  two  reasons.  Firstly, 
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any  nonlinear  disturbance  will  decay  if  forced  at  an  amplitude  that  is  larger  than  its 
natural  nonlinear  saturation  level,  even  if  it  is  strongly  amplified  for  smaller  forcing 
amplitudes.  Secondly,  the  large  spacing  between  the  localized  disturbance  structures 
shown  in  Figure  14,  indicates  that  a  forcing  frequency  of  /  =  110 .8Hz  may  be  far  too 
low  for  allowing  amplified  disturbances.  From  our  experience  with  strong  wall  jets, 
vortical  disturbances  are  strongly  damped  if  the  wave  length  is  much  larger  than  the 
local  length  scale  which,  in  the  present  case,  is  the  wall  jet  thickness  characterized 
either  by  6min  or  Smax.  In  addition,  the  wall  jet  mean  flow  decays  drastically  over  one 
fundamental  disturbance  wave  length  and  the  shape  of  the  weak  wall  jet  profile  with 
its  three  inflection  points,  (one  at  the  wall,  one  in  each  free  shear  layer)  are  essential 
for  the  growth  of  two-dimensional  disturbances.  We  therefore  suspected  that  ampli¬ 
fied  disturbances  exist  for  higher  forcing  frequencies  where  the  wave  length  is  in  the 
order  of  the  wall  jet  thickness.  Therefore,  we  performed  a  parameter  study  to  test 
this  hypothesis. 

From  this  parameter  study,  results  of  four  cases  are  now  presented.  For  all  cases  (case 
1-4),  the  forcing  amplitude  was  chosen  as  Ap  =  0.065  •  U/s  which  is  5%  of  the  steady 
blowing  amplitude  and  ten  times  lower  than  that  of  the  highly  nonlinear  forcing 
discussed  above  (case  1L).  This  should  still  be  small  enough  *o  allow  for  disturbance 
growth  before  nonlinear  saturation  sets  in.  Figure  16  illustrates  the  four  cases  with 
color  contour  plots  of  instantaneous  spanwise  vorticity.  From  case  to  case,  the  forcing 
frequency  is  increased  as  indicated  in  the  figure  captions.  For  case  1,  with  the  lowest 
frequency  (Figure  16a),  no  vortical  structures  are  discernible.  The  forcing  frequency 
matches  that  of  case  1L  (Figure  14),  but  for  the  lower  forcing  amplitude  used  in  case 
1,  the  higher  harmonics  have  a  much  lower  amplitude  and  the  disturbance  is  not 
localized.  As  the  forcing  frequency  is  increased  (Figures  16b-d),  vortical  structures 
become  visible.  They  appear  strongest  for  case  3  in  Figure  16c.  As  expected  from 
the  decay  of  the  mean  flow,  the  disturbances  eventually  decay  for  all  cases. 


a)  case  1 :  f=l  10.8Hz 


b)  case  2:  f=554Hz 


c)  case  3:  f=l  108Hz 


d)  case  4:  f=2216Hz 
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Figure  16:  URANS  calculations  of  pulsed  weak  wall  jet  (medium  grid  resolution). 
Shown  are  color  contours  of  instantaneous  spanwise  vorticity  for  three  different  forcing 
frequencies.  For  all  cases,  steady  wall-tangential  blowing  occurs  with  As  =  1.3  •  Ufs, 
pulsed  blowing  occurs  with  As  =  0.065  •  U fs  while  the  forcing  frequency  varies  as 
indicated. 
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A  more  detailed  analysis  of  the  magnitude  of  the  disturbances  and  their  streamwise 
amplitude  development  is  presented  in  Figure  17.  Shown  for  all  four  cases  are  three 
characteristic  measures  of  the  disturbances  versus  distance  from  the  center  of  the 
blowing  slot.  Figure  17a  depicts  the  fundamental  wave  length  of  the  disturbances  in 
comparison  with  the  wall  jet  thicknesses,  5max  and  5m*n.  Once  traveling  disturbances 
establish  shortly  downstream  of  the  blowing  slot,  their  wave  lengths  stay  fairly  con¬ 
stant  for  all  cases.  The  wave  length  for  case  3  is  about  twice  Smin,  while  for  case  4  it 
approaches  5min.  The  decrease  in  wave  length  from  case  to  case  directly  correlate  with 
the  increase  in  the  forcing  frequency,  which  indicates  that  they  all  travel  at  about  the 
same  phase  speed. 

The  phase  speed  of  the  disturbances  is  shown  in  Figure  17b  in  comparison  with  the 
maximum  and  minimum  velocity  of  the  wall  jet  mean  flow,  Umax  Umin •  For  cascss 
2-4,  the  phase  speed  lies  in  between  Umax  and  Umin  while  for  case  1  it  is  at  or  above 
Umax-  The  disturbance  for  case  1  therefore  decays,  while  growth  is  possible  for  cases 
2-4.  The  phase  speed  of  «  0.8  •  Uj3  also  indicates  that  the  disturbances  are  connected 
with  an  instability  of  the  free  shear  layer,  since  boundary  layer  instabilities  generally 
travel  at  lower  speeds. 

Finally,  Figure  17c  shows  the  amplitude  maximum  of  the  fundamental  disturbance 
on  a  semi-logarithmic  scale.  The  disturbance  for  case  1  indeed  decays  monotonically, 
while  for  cases  2  and  3  initial  growth  is  observed,  although  not  very  strong  and  not  for 
long.  Interestingly,  farther  downstream  the  high  frequency  forcing  cases  2-4  exhibit 
stronger  decay  than  the  low  frequency  case  1. 

In  Figure  18,  the  normalized  shape  of  the  disturbance  amplitudes  are  plotted  for  three 
streamwise  locations.  For  all  disturbances,  the  largest  amplitude  peak  is  in  the  free 
shear  layer  of  the  wall  jet,  clearly  indicating  its  connection  with  that  portion  of  the 
flow  profile. 
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Figure  17:  URANS  calculations  of  pulsed  weak  wall  jet,  cases  1-4  (see  Figure  16a- 
d).  Shown  are  a)  wave  length,  b)  phase  speed,  and  c)  amplitude  of  the  disturbances 
versus  distance  from  the  blowing  slot. 


Figure  18:  URANS  calculations  of  pulsed  weak  wall  jet,  cases  1-4  (see  Figure  16a-d). 
Shown  are  the  Fourier  amplitudes  of  the  streamwise  disturbances  at  three  streamwise 
locations  (a-c)  corresponding  to  those  in  Figure  13a. 
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6  Summary  and  Outlook 

With  the  present  research  we  have  successfully  demonstrated  that  the  problem  of 
separation  control  using  wall  tangential  blowing  can  be  tackled  with  LES,  DNS  and 
a  new  Flow  Simulation  Methodology  (FSM).  Preliminary  calculations  using  a  con¬ 
ventional  LES  have  shown  that  wall  blowing  is  beneficial  for  delaying  boundary  layer 
separation.  The  calculations  also  indicated  that  pulsed  blowing  is  more  effective  than 
steady  blowing.  In  the  LES,  the  pulsed  blowing  lead  to  the  formation  of  2D  co¬ 
herent  structures.  However,  details  of  the  formation  process,  the  evolution  of  these 
structures,  and  the  mechanism  for  delaying  separation  could  not  be  determined  with 
confidence,  since  the  calculations  lacked  the  proper  resolution.  The  appropriate  grid 
resolution  for  the  LES  would  have  been  beyond  the  capacity  of  the  available  compu¬ 
tational  resources. 

A  very  promising  alternative  to  3D  LES  are  our  FSM  calculations  in  the  2D  URANS 
limit.  This  approach  was  proven  to  be  extremely  successful,  when  applied  to  the 
benchmark  case  of  a  strong  turbulent  wall  jet  and  yielded  remarkable  agreement  with 
the  experimental  results.  Since  URANS  calculations  are  two-dimensional,  they  are 
much  less  expensive  than  LES  or  DNS,  and  a  much  highe.  wall-normal  resolution 
can  be  afforded.  For  the  unsteady  weak  wall  jet,  this  enables  to  resolve  the  coherent 
structures  well  and  to  obtain  the  true  shape  of  the  amplitude  distributions  (“eigen¬ 
functions”).  We  have  started  with  parameter  studies  for  a  weak  wall  jet  that  is  close 
to  experiments  by  Wygnanski  and  coworkers.  Our  initial  results  from  these  URANS 
calculations  are  very  encouraging.  We  could  demonstrate  that  vortical  disturbance 
can  be  generated  with  our  approach.  For  smaller  forcing  amplitudes  and  within  a 
certain  frequency  range,  disturbance  growth  could  be  observed,  even  though  only  for 
a  short  distance.  For  large  amplitude  pulsing,  low  frequency  forcing  may  prove  to  be 
the  better  choice,  since  the  disturbances  decay  slower  than  for  higher  frequencies. 
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Many  open  questions  remain,  however.  For  example,  from  our  investigations  it  is 
still  unclear  at  present,  whether  the  weak  wall  jet  allows  for  subharmonic  resonance. 
This  is  a  prominent  mechanism  in  transitional  and  turbulent  strong  wall  jets  and  pro¬ 
duces  far  greater  growth  rates  than  primary  instability  (see  appendix).  Potentially, 
a  subharmonic  cascade  could  be  very  helpful  in  naturally  sustaining  large  coherent 
structures  in  the  flow.  A  second  issue  which  deserves  further  study  is  the  receptivity 
problem.  So  far,  from  the  many  possible  forcing  methods,  only  longitudinal  pulsing 
of  the  slot  velocity  has  been  employed  which  resulted  in  a  shear  layer  mode  of  insta¬ 
bility.  With  a  more  localized  forcing  method  (e.g.  applying  a  local  force  field  inside 
the  flow),  different  instability  modes  could  be  excited  depending  on  the  location  of 
the  forcing.  For  these  detailed  future  investigations,  a  more  realistic  modeling  of 
the  slot  geometry  is  clearly  desirable.  This  is  a  very  difficult  problem,  however,  since 
“real”  blowing  slots  in  experiments  and  in  technical  wall  jet  applications  may  produce 
transitional,  not  fully  turbulent  flow.  In  general,  it  is  extremely  difficult  to  achieve 
a  good  match  with  experiments,  because  of  the  enormous  set  of  parameters  that  can 
influence  the  wall  jet  flow,  and  which  may  not  always  be  readily  available. 
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Abstract.  The  competition  between  2-D  and  3-D  resonances  in  wall  jet  transition  is 
analyzed  using  DNS.  In  2-D  DNS  strong  growth  of  subharmonic  disturbances  is  ob¬ 
served  which  leads  to  a  subharmonic  resonance  cascade.  When  a  periodically  forced 
flow  is  simultaneously  perturbed  by  a  small  amplitude  pulse,  massive  vortex  ejections 
occur.  With  3-D  DNS  it  is  shown  that  additional  3-D  forcing  can  prevent  vortex 
mergings  and  vortex  ejections  by  reducing  the  spanwise  coherence  of  the  vortices,  in 
particular  close  to  the  wall. 

Key  words:  wall  jet,  subharmonic  resonance,  pulse  disturbance,  vortex  ejection 

1  Introduction 

The  transitional  Glauert  wall  jet  [2]  is  studied,  because  this  flow  geometry  is  proto¬ 
typical  for  important  technical  applications,  such  as  film-cooling  and  boundary  layer 
control.  While  the  Glauert  wall  jet  is  geometrically  simpler  than  technical  wall  jet 
flows,  it  exhibits  the  same  main  characteristics,  the  presence  of  a  boundary  layer 
(near- wall  region)  in  close  proximity  to  a  free  shear  layer  (outer  region).  Due  to  the 
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inflection  point  in  the  velocity  profile,  the  outer  region  of  the  wall  jet  is  susceptible  to 
an  inviscid  Kelvin-Helmholtz  type  instability  which  may  give  rise  to  two-dimensional 
(2-D)  subharmonic  resonances.  Due  to  the  presence  of  the  wall,  the  near-wall  region 
is  governed  by  viscous  instability.  Thus,  growth  of  Tollmien-Schlichting  type  waves 
can  occur  allowing  for  three-dimensional  (3-D)  secondary  instabilities  (fundamental 
and  subharmonic  resonances).  Therefore,  in  a  wall  jet  2-D  and  3-D  resonance  mech¬ 
anisms  can  exist  simultaneously  and  may  compete  with  each  other.  In  fact,  various 
transition  scenarios  have  been  observed  in  experiments  of  forced  transitional  wall  jets. 
For  example,  in  experiments  by  Bajura  and  Catalano  [1],  2-D  subharmonic  resonance 
plays  a  major  role  during  transition,  whereas  in  experiments  by  Shih  and  Gogineni  [4], 
three-dimensionality  sets  in  immediately  downstream  of  the  primary  disturbances. 

In  previous  work,  using  direct  numerical  simulations  (DNS),  we  have  investigated  2-D 
secondary  instability  [5]  and  3-D  secondary  instabilities  [6]  separately.  For  the  inves¬ 
tigated  cases,  we  concluded  that  3-D  subharmonic  resonance  is  at  least  as  strong  as 
the  2-D  subharmonic  resonance.  With  the  present  study  (also  using  DNS),  transition 
scenarios  are  investigated  where  both  2-D  and  3-D  resonances  occur  simultaneously. 
Through  selective  forcing  the  conditions  are  explored  under  which  2-D  resonances  or 
3-D  resonances  dominate  the  transition  process. 


2  Computational  Approach 

For  the  DNS,  a  3-D  incompressible  Navier-Stokes  code  for  DNS  of  boundary  layer 
transition  [3]  has  been  adapted  to  the  wall  jet  geometry.  In  the  DNS  the  incompress¬ 
ible  vorticity  transport  equations  (VTE)  and  velocity  Poisson  equations  are  solved 
in  total  flow  formulation  inside  the  computational  domain  shown  in  Fig.  1.  The  nu¬ 
merical  scheme  employs  a  fourth-order  accurate  Runge-Kutta  method  for  the  time 
integration  and  fourth-order  compact  differences  in  streamwise  and  wall-normal  di¬ 
rections.  For  the  wall-normal  direction  a  variable  grid  is  used  with  points  clustered 
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Qlausit  hase  flow: 

kinematic  viscosity: 
\=1.5xlO~sm2/sec 

flux  of  exterior 
momentum  flux: 

F=I.68x  IO~5ms/sec 3 

local  inflow 
Reynolds  number: 

ReSos=280  at  x=5 

Figure  1:  Computational  domain  with  Glauert  base  flow  and  blowing/suction  slot  for 
disturbance  generation 

near  the  wall.  In  3-D  DNS  the  spanwise  direction  is  treated  pseudo-spectrally.  The 
boundary  conditions  include  Dirichlet  conditions  at  the  inflow,  zero  disturbances  at 
the  free  stream  boundary,  and  a  buffer  domain  near  the  outflow  boundary  for  relam- 
inarizing  the  flow  [3].  At  the  wall,  the  velocity  is  set  to  zero  except  over  the  blowing 
and  suction  slot.  In  a  DNS~precursor  calculation,  a  2-D  base  flow  is  computed  (Fig.  1) 
which  closely  matches  the  simi-  larity  solution  that  was  used  in  [5]  as  a  base  flow.  To 
ensure  that  the  DNS  base  flow  is  independent  of  the  computational  domain  height, 
special  care  must  be  taken  to  recover  the  induced  flow  field  of  the  wall  jet  far  away 
from  the  wall,  where  the  similarity  solution  is  invalid.  This  is  achieved  by  imposing 
flow  profiles  at  the  inflow  that  combine  the  similarity  solution  near  the  wall  with  an 
asymptotically  matched  potential  flow  solution  away  from  the  wall. 

2.1  Forcing  Method 

For  the  transition  simulations,  the  base  flow  is  forced  by  simultaneous  blowing  and 
suction  through  a  slot  in  the  wall  close  to  the  inflow  boundary.  This  technique  is  very 
efficient  in  producing  vortical  disturbances  while  minimizing  acoustic  disturbances. 
The  slot  velocity  is  computed  as 

K  N 

vs{x,  y  =  0,z,  t )  =  FVt{x)  ■  J2  Y  A(n*)  ■  cos(2tt fnt)  •  cos(t kz) 

k=0  n=0 


(14) 
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Figure  2:  Forcing  method  (a)  streamwise  variation  FVs(x )  (b)  mode  diagram  for 
forcing  amplitudes  ^4(n,Jt)  of  periodic  forcing  (c)  amplitude  Ap(t )  of  pulse  disturbance 
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1 
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a) 0  1  n  3  4  b)  0  1  S  3  4  C) 

Figure  3:  Mode  diagrams  for  forcing  amplitudes  (a)  Case  1  (b)  Case  2  (c)  Case  3 

where  A(n<k)  is  the  forcing  amplitude  of  mode  (n,  k )  in  [%],  /„  =  n-  f  is  the  frequency 
in  [Hz],  7 k  =  k  ■  2-k/\z  is  the  spanwise  wave  number  in  [1/3. 2mm].  The  forcing 
method  is  illustrated  in  Fig.  2.  At  first,  the  flow  is  forced  periodically  (Fig.  2b)  until 
it  converges  to  a  time-periodic  state.  Then,  at  time  tp,  a  small  pulse  disturbance  (Fig. 
2c)  is  also  introduced  in  mode  (4,0).  In  the  present  paper,  three  cases  with  different 
periodic  forcing  are  discussed  (Fig.  3).  The  pulse  disturbance  (Fig.  2c)  is  identical 
for  all  cases. 


2.2  Computational  Parameters 

The  flow  quantities  are  non-dimensionalized  with  a  reference  length,  L  =  3.2mm,  and 
the  base  flow  velocity  at  the  slot,  U  =  3m/ sec.  For  the  discretization  in  streamwise 
direction,  a  constant  stepsize,  Aa;  =  0.125L,  is  employed  for  all  cases.  The  domain 
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length  varies  however:  For  Case  1,  1201  points  (+  800  points  in  buffer  domain)  are 
used.  For  Case  2,  1201  (+  400)  points  are  used,  and  for  Case  3,  801  (+  200)  points. 
In  wall-normal  direction,  an  exponentially  stretched  grid  with  180  points  is  employed 
where  the  smallest  stepsize,  A ym)n  =  0.02L,  is  at  the  wall,  and  the  largest  stepsize, 
A ymax  =  1.2 L,  is  at  the  upper  boundary.  The  spanwise  direction  is  composed  of  11 
symmetric  spectral  modes  with  the  largest  spanwise  wave  number  yi  =  0.1  ■!hi/L. 


2.3  Validation  of  Computations 

For  validation  of  the  DNS  code  see  Meitz  [3].  For  the  DNS  of  wall  jets,  convergence 
investigations  regarding  the  effect  of  stepsizes  have  ensured  that  the  the  influence 
of  the  discretization  error  on  the  flow  solution  is  kept  small.  Especially  for  2-D 
simulations,  a  long  buffer  domain  is  crucial  for  preventing  upstream  feed  from  the 
outflow  boundary.  Even  small  feed  back  can  excite  strongly  amplified  low  frequency 
disturbances  which  can  prevent  the  flow  from  reaching  a  time-periodic  state. 
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3  Results 

We  have  conducted  numerical  experiments  to  uncover  the  relative  importance  of  2-D 
and  3-D  secondary  instabilities.  Three  case  studies  are  presented  here. 

3.1  2-D  Flow  -  Case  1 

In  Case  1,  the  flow  is  forced  with  a  fundamental  and  two  subharmonic  frequencies,  all 
2-D  (Fig.  3a).  The  resulting  time-periodic  flow  exhibits  two  subsequent  vortex  merg¬ 
ings  (Fig.  4a).  A  Fourier  decomposition  of  the  time-periodic  flow  (Fig.  4b)  reveals 
that,  through  resonance  with  the  large  fundamental  disturbances,  the  small  subhar¬ 
monic  disturbances  are  strongly  amplified  and  successively  surpass  the  amplitude  of 
the  fundamental.  This  subharmonic  cascade  is,  however,  extremely  sensitive  to  small 
non-periodic  disturbances. 

When  at  time  tp  the  wall  jet  is  additionally  perturbed  with  a  small  pulse  (Fig.  2c), 
it  responds  by  ejecting  dipolar  vortices  away  from  the  wall  (Fig.  5a)  and,  in  fact 
completely  lifts  off  from  the  wall.  As  illustrated  in  the  x  -  t  diagram  in  Fig.  5b, 
the  pulse  disturbance  generates  a  downstream  traveling  wave  packet  which  grows  in 
amplitude  by  two  orders  of  magnitude  and  exhibits  several  doublings  of  the  dominant 
wave  length.  In  contrast,  without  the  presence  of  large,  oeriodic  disturbances  (not 
shown),  the  wave  packet  barely  grows  and  experiences  only  a  gradual  shift  to  smaller 
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Figure  4:  Time-periodic  flow  -  Case  1  (a)  grey  scales  of  instantaneous  wz-vorticity 
(b)  Fourier  amplitudes  of  u;z-vorticity  at  y  =  0  for  three  modes 
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Figure  5:  Flow  response  to  small  pulse  disturbance  —  Case  1  (a)  grey  scales  of  u)z- 
vorticity  (b)  x-t  diagram  of  wz-vorticity  at  y  =  0  (flow  with  pulse  -  without  pulse) 

wave  numbers.  This  indicates  that  the  subharmonic  components  in  the  wave  packet 
are  strongly  amplified  through  resonance  with  the  large  amplitude  periodic  distur¬ 
bances.  As  the  wave  packet  reaches  the  amplitude  level  of  the  periodic  disturbances 
it  disrupts  the  balance  of  the  double  vortex  row  and  causes  the  ejection  of  vortex 
pairs. 

3.2  3-D  Flow  —  Cases  2  and  3 

For  Case  2,  in  addition  to  the  disturbances  of  Case  1,  3-D  disturbances  with  the  same 
amplitudes  as  the  2-D  subharmonic  disturbances  are  introduced  (Fig.  3b).  Now,  the 
flow  becomes  highly  three-dimensional  downstream  of  x  =  70.  As  illustrated  in  Fig. 
6a,  the  same  pulse  as  in  Case  1  causes  no  vortex  ejections,  but  merely  produces  a 
number  of  small  vortices  that  linger  in  the  outer  region  of  the  wall  jet  before  they  are 
convected  out  of  the  computational  domain.  How  differently  the  wave  packet  now 
develops  compared  to  Case  1  becomes  evident  from  the  x  —  t  diagrams  for  the  two 
cases  (Figs.  6b  and  5b).  Initially,  both  cases  are  virtually  identical.  But  downstream 
of  x  =  70,  where  in  Case  1  the  vortex  ejection  occurs,  in  Case  2  the  amplitude  of  the 
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Figure  6:  Flow  response  to  small  pulse  -  Case  2  (a)  grey  scales  of  u;z-vorticity  at 
z  =  0  (b)  x-t  diagram  of  u;z-vorticity  for  y  =  0,  k  =  0  (with  pulse  -  without  pulse) 


Figure  7:  Flow  response  to  small  pulse  -  Case  3  (a)  grey  scales  of  wz-vorticity  at 
z  =  0  (b)  x-t  diagram  of  uz-vorticity  for  y  =  0,  k  =  0  (with  pulse  -  without  pulse) 


spectral  component  k  =  0  of  the  wave  packet  decreases  markedly.  The  wave  packet 
looses  its  spanwise  coherence  along  with  the  periodic  2-D  disturbances. 

Only  when,  as  in  Case  3,  the  amplitude  of  both  the  2-D  fundamental  and  subharmonic 
forcing  is  significantly  increased  (Fig.  3c),  vortex  merging  and  vortex  ejection  occur 
prior  to  development  of  the  three-dimensional  stage  (Fig.  7a).  Now  the  wave  packet 
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Figure  8:  Flow  response  to  small  pulse  disturbance.  Shown  are  iso-surface  plots  of 
u;z-vorticity  at  time  t  =  tp  +  30  (a)  Case  2  (b)  Case  3. 

reaches  a  large  amplitude  level  and  exhibits  a  doublings  of  the  dominant  wave  length 
(Fig.  7b). 

The  difference  in  the  transition  process  between  Case  2  and  Case  3  and  the  impact  of 
the  small  pulse  in  each  case  is  nicely  illustrated  by  the  3-D  snapshots  in  Fig.  8.  While 
in  Fig.  8a  the  flow  becomes  highly  three-dimensional  immediately  downstream  of  the 
fundamental  disturbance,  in  Fig.  8b  vortex  mergings  and  the  ejection  of  a  dipolar 
vortex  precede  the  3-D  stages. 

4  Conclusion 

In  our  2-D  DNS  of  a  forced  wall  jet,  we  found  that  subharmonic  resonance  is  the 
dominant  mechanism.  It  strongly  amplifies  small  (subharmonic)  disturbances  and 
thus  renders  the  flow  extremely  sensitive  to  small  amplitude  pulse  (or  random)  dis¬ 
turbances  causing  vortex  ejections  and  the  possible  disintegration  of  the  wall  jet.  In 
our  3-D  DNS,  the  2-D  subharmonic  resonance  may  be  surpassed  by  3-D  resonances 
which,  as  a  consequence,  inhibit  vortex  mergings  and  vortex  ejections.  Our  findings 
may  explain  experimental  observations,  where  vortex  mergings  and  vortex  ejections 
occur  for  low  Reynolds  numbers  and  low  turbulence  levels  [1],  but  not  for  higher 
Reynolds  numbers  and  higher  turbulence  levels  [4]. 
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Abstract 

The  effect  of  time  periodic  forcing  of  turbulent  wall 
jets  is  investigated  numerically  using  our  new  Flow 
Simulation  Methodology  (FSM).  Previous  experi¬ 
mental  research  has  shown  that  periodic  excitation 
of  wall  jets  increases  their  effectiveness  for  separa- 
tion  control.  To  investigate  the  effect  of  forcing 
on  the  large  coherent  structures,  a  steady  mean 
how  is  computed  first  using  a  Reynolds  Averaged 
Xavior-Stokos  (RANS)  approach.  The  response  of 
this  mean  how  to  periodic  forcing  is  then  investi¬ 
gated  using  FSM  in  the  limit  of  unsteady  RANS. 
Our  results  indicate  that  with  periodic  forcing,  sig¬ 
nificant  changes  in  the  mean  Reynolds  shear  stress 
distribution  and  also  the  skin  friction  are  observed. 
The  unsteady  computational  results  are  scrutinized 
to  investigate  how  the  coherent  structures  interact 
with  the  turbulent  mean  how  and  what  role  they 
plav  m  the  change  of  the  mean  how  profiles. 

1  Introduction 

A  wall  jet  develops  when  a  fluid  is  blown  tangen¬ 
tially  along  a  surface.  Wall  jets  have  many  impor¬ 
tant  technical  applications  such  as  boundary  layer 
control  and  lift  enhancement  on  airfoils,  making 
use  of  their  tendency  to  adhere  to  curved  surfaces 
i Coanda  effect).  In  other  applications,  such  as  the 
cooling  of  turbine  blades  and  combustion  chamber 
walls,  a  wall  jet  is  used  to  cool  and  protect  surfaces 
from  hot  and/or  corrosive  fluids. 

Irrespective  of  its  application,  a  wall  jet  acts  as  a 
source  of  momentum.  The  jet  may  be  issued  into  a 
eohowing.  quiescent,  or  counterhowing  surrounding. 
The  external  stream  may  be  accelerating  or  deceler- 
ating  such  as  in  airfoil  or  turbine  blade  applications. 
In  addition,  the  external  how  can  be  laminar  or  tur- 
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bulent.  The  surface  itself  might  be  smooth,  rough, 
or  even  highly  irregular,  such  as  in  electronic  cool¬ 
ing  applications.  Despite  its  applications  in  many 
technical  systems,  the  wall  jet,  especially  the  role  of 
large  coherent  structures  in  the  turbulent  wall  jet. 
is  not  at  all  understood. 

Experimental  work  by  Wygnanski  and  cowork¬ 
ers  [1]  has  demonstrated  the  effectiveness  of  using 
pulsed  wall  jets  for  delaying  separation.  Using  flow 
visualization,  Wygnanski  was  able  to  show  large  co¬ 
herent  structures  in  the  turbulent  flow  on  an  air¬ 
foil.  These  structures  profoundly  affect  the  mean 
flow,  significantly  delaying  separation.  More  im¬ 
portantly,  deliberately  introducing  these  structures 
using  periodic  forcing  results  in  predictable  changes 
of  the  mean  flow. 

The  geometry  used  in  the  experiments  is  highly 
complex,  combining  the  effects  of  wall  curvature, 
adverse  pressure  gradient,  a  coflowing  external 
stream,  and  unsteady  momentum  addition  by  the 
pulsed  wall  jet.  To  understand  the  most  dominant 
physical  mechanisms  in  this  flow,  it  is  paramount  to 
understand  the  influence  of  every  component  of  the 
flow  in  the  experimental  investigation.  As  a  first 
step  in  a  series  of  numerical  investigations,  the  ef¬ 
fect  of  periodic  forcing  on  the  turbulent  mean  flow 
of  the  strong  turbulent  wall  jet  on  a  flat  surface  is 
studied  in  this  paper. 

A  large  number  of  experimental  investigations  of 
the  mean  flow  of  turbulent  wall  jets  has  been  re¬ 
ported  in  the  literature  (see  review  articles  [2],  [3], 
[4]).  Although  the  measurements  are  very  consis¬ 
tent  in  a  large  part  of  the  mean  profiles,  the  velocity 
in  the  near  wall  region  is  extraordinarily  difficult  to 
measure  and  therefore  skin  friction  data  vary  con¬ 
siderably.  Furthermore,  since  these  experimental 
investigations  focused  on  the  mean  flow  quantities 
only,  the  question  of  the  existence  and  influence  of 
large  coherent  structures  in  the  outer  shear  layer  on 
the  mean  flow  profiles  was  not  addressed. 

In  this  paper,  the  role  of  these  large  coher¬ 
ent  structures  is  investigated  numerically  using  our 
new  Flow  Simulation  Methodology  (FSM).  This  ap- 
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proach  has  been  developed  in  collaboration  with  C. 
Speziale  [5]  because  of  the  limitations  in  traditional 
computational  strategies  such  as  Direct  Numeri¬ 
cal  Simulation  (DNS),  Reynolds-averaged  Navier- 
Stokes  (RANS)  computation,  and  Large  Eddy  Sim¬ 
ulation  (LES).  DNS  is  limited  to  relatively  low 
Reynolds  number  flows  and  is,  in  the  foresee¬ 
able  future,  impractical  for  high  Reynolds  num¬ 
ber  flows  which  are  of  greatest  engineering  inter¬ 
est.  This  leaves  Reynolds-averaged  Navier-Stokes 
(RANS)  calculations  and  Large-Eddy  Simulation 
(LES)  as  the  only  realistic  alternatives  to  compute 
turbulent  flows  [5].  However,  traditional  RANS 
methods  are  limited  to  time-steady  solutions,  an 
assumption  unsuitable  for  complex  turbulent  flows 
where  unsteady  large  structures  play  a  dominant 
role  in  the  transfer  of  heat  and  momentum  between 
different  regions  of  the  flow.  From  a  theoretical 
standpoint,  LES  has  long  been  viewed  as  an  ideal 
method  for  the  time-dependent,  spatially  varying 
solution  of  complex  turbulent  flows.  Unfortunately, 
the  promise  of  LES  to  provide  reasonable  solutions 
for  larger  Reynolds  numbers  remains  largely  unful¬ 
filled. 

FSM  is  aimed  at  combining  the  advantages  of 
RANS  and  LES  in  a  consistent  manner.  In  fact, 
FSM  maintains  the  general  form  of  LES,  i.e.,  the 
large  scales  are  computed  and  the  small  scales  are 
modelled.  The  main  feature  of  our  FSM  is  that 
the  subgrid-scale  stress  model  of  LES  is  replaced 
by  t  he  product  of  a  contribution  function  and  a 
state-of-the-art  Reynolds  stress  model.  The  con¬ 
tribution  function  is  designed  such  that,  depend¬ 
ing  on  the*  numerical  resolution  (as  compared  to  a 
turbulent  length  scale  like  the  Kolmogorov  length 
scale),  the  simulations  approach  a  DNS  in  the  limit 
of  fine  resolution  (or  low  Reynolds  number),  or  an 
unsteady  RANS  in  the  limit  of  coarse  resolution  (or 
high  Reynolds  number). 

Using  FSM  in  the  unsteady  RANS  limit,  the  role 
of  the  large  coherent  structures  in  the  turbulent 
wall  jet  is  investigated  in  this  paper.  Experimen¬ 
tal  investigations  ([6],  [7])  have  shown  that  the  skin 
friction  is  markedly  reduced  when  large  amplitude 
disturbances  are  introduced  into  the  turbulent  wall 
jet.  Furthermore,  the  experiments  by  Schober  [8] 
show  that  the  resulting  coherent  structures  are  pre¬ 
dominantly  two-dimensional.  If  these  structures  are 
in  fact  two-dimensional,  a  two-dimensional  simula¬ 
tion  in  the  FSM  limit  of  unsteady  RANS  should  be 
able  to  capture  the  most  important  effects  of  the 
coherent  motion  on  the  turbulent  mean  flow. 


Figure  1:  Computational  domain  with  schematic  of 
volume  forcing  for  disturbance  generation 

2  Numerical  Method 

2.1  Governing  equations 

The  governing  equations  are  the  two  dimensional, 
incompressible,  unsteady  Navier-Stokes  equations 
in  vorticity- velocity  formulation.  The  velocity  com¬ 
ponents  in  the  streamwise  (:r)  and  wall  normal  (y) 
directions  are  u  and  v,  respectively  (figure  1). 

In  these  equations,  the  velocities  are  normalized 
by  the  jet  exit  velocity  Uj,  the  spatial  variables  x,  y 
by  the  nozzle  width  6,  and  the  time  by  b/Uj.  Taking 
the  curl  of  the  momentum  equations,  which  elimi¬ 
nates  the  pressure  gradient  terms,  yields  the  trans¬ 
port  equations  for  the  vorticity, 

^  +  Vx(3xr-V-?)  =  d-V23.  (1) 

dt  V  /  Rej 

where  the  overbar  denotes  an  averaging  operation 
(temporal  or  spatial),  f  =  [r^]  is  the  turbulent 
stress  tensor.  In  the  two-dimensional  case,  this  set 
of  equations  reduces  to  one  equation  fov  the  span- 
wise  vorticity  component  u)z.  The  velocity  compo¬ 
nents  are  computed  from  a  set  of  Poisson  equations 
(see  [9]  for  more  details). 

In  order  to  close  equation  (1),  the  turbulent  stress 
tensor  components  rXJ  need  to  be  modelled.  The 
model  used  in  this  investigation  is  our  new  FSM. 
Following  Speziale  [5],  the  turbulent  stress  tensor 
components  are  written  as 

TlJ  =  /(A/L,)r«  ,  (2) 

where  r-j  is  the  Reynolds  stress  tensor  and  /( A/L*) 
is  the  contribution  function.  In  the  contribution 
function,  A  is  the  effective  computational  grid  size 
and  Lk  =  [v2  A) 1/4  is  the  Kolmogorov  length  scale. 
Note  that  the  contribution  function  depends  on  the 
ratio  A/L*),  i.e.,  the  physical  resolution  of  the  com¬ 
putation  compared  to  a  length  scale  of  turbulence 
in  the  flow.  The  Reynolds  stress  r in  equation  (2) 
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is  computed  using  the  anisotropic  Algebraic  Stress  Finally,  the  constants  needed  for  the  solution  of 

Model  (ASM)  of  Gatski  and  Speziale  [10],  equations  (7)  through  (9)  are, 


-mt) 


k2^ 


or 


+a2^{SikWkj  +  SjkWkl) 

~ *-*3^7  ^SikSkj  ~  klS kl^ij^ 


(3) 


The  strain  rate  tensor  and  vorticity  tensor  are 
given  by 

_UduL,dul\w  i  fdui_dul\ 

b,J  -  2\dxj  +  dxij’  ij  2  \dxj  dxij- 

(4) 

The  function  /(r;,  £ ),  in  conjunction  with  k2/e, 
assumes  the  role  of  a  nonlinear  eddy  viscosity, 


/(»/•  0  = 


3(1  +  rj2) 


3  -  2 t]2  +  6£2  ~  3  +  T}2  +  6f?2£2  +  6f2  ’ 

(5) 

where  the  expression  on  the  right  hand  side  of  equa¬ 
tion  5  is  a  regularized  form  of  f(i], f),  which  is  used 
to  avoid  a  possible  division  by  zero  in  the  original 
expression  [5]. 

1 1  and  £  are  invariants  of  the  irrotational  strain 
rate  tensor  and  the  vorticity  tensor,  respectively, 


1  «3  k  ,1/2  c  _  02  k  /TTT  TJT  ,1/2 

t/-9n  (S'jS.j)  ,£  -  ai  £  (W  ijW'y)  • 

(6) 

The  values  of  model  constants  in  equations  (3) 
and  (G)  are  Oj  —  0.227,  Q2  =  0.0423,  and  03  = 
0.039G. 

The  turbulent  kinetic  energy  k  and  the  turbulent 
dissipation  rate  e  in  equation  (3)  are  computed  from 
the  standard  transport  equations  [11], 


Dk  dux  d 

Dt  TlJ  dxj  t  +  dxj 


i/'Y  1  \  dk 

Gk  +  Re )  dxj  J 


[7) 


D(  -  r  eT  —  r  f  - 

Di~  C<lk,JdxJ  C(2f<2k 


+- 


l/T  1 

—  + 


dxj  [  \ae  Re  J  dx3 


8e 


•(B) 


The  eddy-viscosity  ist  and  the  wall  damping  func¬ 
tion  }t2  are  defined  as 


vt  —  C^k2/e ,  /t2  =  l-e-s'fle^/1°.  (9) 


C„  =  0.09 ,  Cei  =  1.44,  C(  0  =  1.83, 

=  1 ,  a<  —  1.3,  -4+  =  25. 

2.2  Numerical  scheme 

A  detailed  description  of  the  numerical  schemes 
used  for  the  solution  of  equation  (1)  and  the  Pois¬ 
son  equations  for  the  velocities  is  given  in  [9].  For 
the  time  integration,  a  fourth-order  accurate  ex¬ 
plicit  Runge-Kutta  method  is  used.  The  deriva¬ 
tives  in  the  streamwise  direction  are  approximated 
by  fourth-order  compact  upwind  and  downwind  bi¬ 
ased  differences  at  consecutive  time  steps  and  wall- 
normal  derivatives  are  computed  with  fourth-order 
accurate  compact  differences.  Lastly,  a  feist  Fourier 
solver  is  implemented  for  the  solution  of  the  velocity 
Poisson  equations. 

In  addition  to  the  solution  of  the  vorticity  trans¬ 
port  equation  and  the  velocity  Poisson  equations, 
FSM  requires  the  simultaneous  integration  of  the  k 
and  e  equations.  For  these  equations,  a  second  or¬ 
der  ADI  method  is  used  for  the  time  integration  and 
second-order  accurate  difference  stencils  are  em¬ 
ployed  for  the  streamwise  and  wall  normal  deriva¬ 
tives.  The  fourth-order  accurate  difference  method 
used  in  the  discretization  of  the  vorticity  transport 
equation  is  employed  for  the  spatial  derivatives  of 
the  turbulent  stress  tensor  components  (see  equa¬ 
tion  (1)). 

2.3  Initial  and  Boundary  conditions 

The  governing  equations  are  solved  in  the  compu¬ 
tational  domain  shown  in  Figure  1.  At  the  inflow 
boundary,  Dirichlet  conditions  are  imposed  for  the 
vorticity,  the  velocity  components,  the  turbulent  ki¬ 
netic  energy  and  the  turbulent  dissipation  rate.  For 
the  computations  presented  here,  the  inflow  condi¬ 
tions  are  obtained  from  a  solution  of  the  boundary 
layer  equations  for  the  turbulent  wall  jet.  At  the 
wall  y  =  0,  the  no-slip,  no-penetration  conditions 
are  imposed  on  the  velocity.  The  turbulent  kinetic 
energy  is  set  to  zero,  and  a  Neumann  condition  is 
used  for  the  turbulent  dissipation  rate.  At  the  free 
stream  boundary  y  =  ymQx ,  a  decay  condition  is 
imposed  on  v,  and  the  wall  normal  derivative  is  set 
to  zero  for  all  other  quantities.  This  requires  the 
domain  height  to  be  large  enough  so  that  all  distur¬ 
bances  have  decayed  sufficiently  and  the  imposed 
boundary  condition  does  not  influence  the  solution. 
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Near  the  outflow  boundary,  a  damping  region  sim¬ 
ilar  to  the  one  proposed  by  [12]  is  used.  At  the 
outflow  boundary  itself,  all  second  derivatives  are 
set  to  zero. 

As  initial  condition,  the  vorticity,  both  velocity 
components,  the  turbulent  kinetic  energy  and  the 
turbulent  dissipation  rate  are  prescribed.  These 
quantities  are  obtained  from  the  solution  of  the 
boundary  layer  equations  for  the  turbulent  wall  jet. 

Disturbances  are  introduced  into  the  flow  using 
a  volume  forcing  technique  [13].  With  this  method, 
predominantly  vortical  disturbances  can  be  gener¬ 
ated  anywhere  in  the  flow  field.  In  the  cases  pre¬ 
sented  in  this  paper,  the  volume  force  is  applied 
near  the  inflow  boundary  and  extends  through  most 
of  the  wall  jet  in  the  wall  normal  direction. 

3  Results 

3.1  Validation 

Two-dimensional  computations  of  the  mean  flow  of 
the  turbulent  wall  jet  are  used  to  validate  our  new 
FSM  in  the  limit  of  traditional  (steady)  RANS.  For 
the  results  shown  here,  the  jet  exit  velocity  is  Uj  = 
30  ms"1  and  the  nozzle  height  is  b  =  5  mm,  which 
yields  a  jet  exit  Reynolds  number  Re3  —  10000. 

The  computational  domain  size  is  40  <  x  <  240, 
Ihnaj  =  GO.  401  points  are  used  in  the  streamwise 
x-dirertion.  yielding  dx  =  0.5.  With  a  timestep  of 
dt  -  0.022.  the  CFL  number  becomes  CFL  =  0.044. 
The  wall  normal  direction  is  discretized  using  200 
points  Grid  points  are  clustered  near  the  wall 
which  results  in  dy  =  0.00234  at  the  wall.  A  grid 
refinement  study  with  the  boundary  layer  equation 
solver  showed  that  this  grid  spacing  gives  grid  in¬ 
dependent  results. 

The  computational  results  are  compared  with  ex¬ 
perimental  data  of  Eriksson  [14].  The  experiments 
were  conducted  in  a  water  tank.  In  the  experiments, 
the  jet  exited  at  a  velocity  of  Uj  =  1  ms”1  through 
a  nozzle  b  —  10  mm  high,  which  yields  a  jet  exit 
Reynolds  number  of  Re3  =  10000.  For  clarity,  only 
the  downstream  location  of  x/b  =  140  is  shown  in 
the  following  plots.  These  profiles  are  representa¬ 
tive  for  a  large  section  of  the  integration  domain, 
excluding  the  regions  close  to  the  inflow  and  out¬ 
flow  boundaries. 

The  computed  mean  profile  for  the  u-velocity  is 
shown  in  figure  2  in  outer  coordinates.  Overall,  very 
good  agreement  with  experimental  data  is  achieved. 
The  velocity  maximum  is  located  at  y/y i/2  =  0.13. 
slightly  closer  to  the  wall  than  the  experimental 


Figure  2:  Streamwise  velocity  in  outer  coordinates, 
o  experimental  results  at  x/b  =  70  [14] 


Figure  3:  Streamwise  velocity  in  wall  coordinates. 

- computational  results,  o  experimental  results 

at  x/b  =  70,  •  at  x/b  ~  150  [14].  •  •  •  logarithmic 
law  [17] 

findings.  When  the  u-velocitv  is  plotted  in  inner  co- 
ordinates  (iC  =  u/uT,  =  yur/v.  uT  =  \frw/p . 
figure  3),  it  is  seen  that  the  viscous  sublayer  is 
resolved  with  about  5  points  within  y+  <  5  and 
the  results  show  excellent  agreement  with  the  theo¬ 
retical  curve  u+  =  y+ .  In  the  logarithmic  layer, 
u+  =  A\ogy+  *F  B,  the  computed  maximum  u- 
velocity  is  slightly  lower  than  the  measured  value 
and  the  theoretical  prediction,  a  tendency  which 
was  observed  for  other  k  —  e  models  (see  [15]  and 
[16]). 

Figure  4  shows  the  wall  jet  half  width  as  a  func¬ 
tion  of  downstream  distance.  The  downstream  lo¬ 
cation  is  corrected  for  the  location  of  the  “virtual 


4 

American  Insitute  of  Aeronautics  and  Astronautics 


J.  Seidel,  H.  Fasel 


AIAA  2000-2317 


nozzlev .  which  is  determined  such  that  the  linear 
upstream  extrapolation  of  ( Uj/Um )2  is  equal  to  one 
for  {x  -  xo)/b  =  0.  In  the  case  presented  here, 
x0/b  ^  23.  The  data  are  scaled  using  the  nor¬ 
malization  proposed  in  [18],  where  the  length  scale 
is  v2 /M  and  M  =  U2b  is  the  jet  momentum  at 
the  nozzle  exit.  Excellent  agreement  with  the  ex¬ 
perimental  results  is  achieved.  The  spreading  rate 
dyi/o/dx  =  0.073  matches  the  value  given  in  [2]. 

The  decay  of  the  maximum  jet  velocity  Um  is 
shown  in  figure  5,  again  normalized  with  the  veloc¬ 
ity  scale  U  =  M/v.  The  results  are  in  very  good 
agreement,  although  the  computational  results  are 
slightly  higher  than  the  curve  fit  of  the  experimen¬ 
tal  data.  A  similar  trend  is  observed  in  Schober  [8] 
when  his  measurements  are  compared  with  the  data 
by  Wygnanski  et  al  [18]. 

In  figure  6  the  skin  friction  coefficient  is  plotted 
;t>  a  function  of  downstream  distance.  The  compu¬ 
tational  results  indicate  a  slightly  lower  decay  rate 
than  the  curve  fit  in  [18].  This  yields  skin  friction 
value's  which  are  approximately  50%  too  high  to¬ 
wards  the  downstream  end  of  the  computational 
domain.  Schober  [8]  pointed  out  that  his  results 
are  also  in  much  better  agreement  with  early  re- 
Milt >  by  Bradshaw  k  Gee  [19],  who  found  a  higher 
skin  friction  than  [18].  Unfortunately,  the  results 
in  [19]  cannot  be  converted  to  the  scaling  of  fig¬ 
ure  G.  Therefore,  figure  7  shows  the  local  skin  fric¬ 
tion  coefficient  cf  =  ru,/ {1/2 pU^)  as  a  function  of 
local  maximum  Reynolds  number  Rem  =  Umym/v- 
In  these  coordinates,  the  computational  results  are 
onlv  about  20%  higher  than  the  curve  fits  reported 
m  [19  and  [20],  and  the  streamwise  development  is 
in  good  agreement  with  their  data.  A  similar  de¬ 
viation  was  reported  for  the  k  —  i  model  used  in 
[21[ 

Results  for  the  Reynolds  shear  stress  R\2  are 
shown  in  both  inner  and  outer  coordinates  in  fig- 
u i  os  8  and  9.  Excellent  agreement  between  experi¬ 
ments  and  computations  is  achieved.  Near  the  wall, 
#//*/,  t  <  0.4,  the  experimental  results  are  slightly 
underpredicted.  This  is  observed  more  clearly  in 
figure  9,  where  R\2  is  plotted  in  wall  coordinates. 
Comparing  the  experimental  results  in  figure  9  and 
figure  3  shows  the  well  known  nonequilibrium  na¬ 
ture  of  the  wall  jet,  indicated  by  the  displacement 
of  the  location  of  zero  Reynolds  shear  stress  from 
the  location  of  maximum  velocity.  An  equilib¬ 
rium  turbulence  model  like  the  ASM  model  used 
in  this  investigation  cannot  predict  this  effect  since 
R\  >  oc  S ij.  This  could  explain  the  deviation  of 
the  computed  R\2  profiles  from  the  experimental 


\ 

\ 


1  io  "  .  * .  100  1000 

(x-x^M/v*  (xl  O'*) 

Figure  4:  Normalized  jet  half  width  as  a  function 

of  normalized  downstream  distance. - curve  fit 

by  [18],  •  •  •  computational  results 


io  100  1000 

{x-xJM/v*  (xl  0  *) 

Figure  5:  Normalized  local  maximum  velocity  as  a 

function  of  normalized  downstream  distance. - 

curve  fit  by  [18],  •  •  •  computational  results 

results.  In  addition,  this  could  also  affect  the  lo¬ 
cation  of  the  velocity  maximum,  which  is  slightly 
closer  to  the  wall  {ym/y  1/2  —  0-13)  than  the  aver¬ 
age  of  experimental  findings,  ym/y\/2  %  0.15  (see 
inset  in  figure  2). 

Figures  10  and  11  show  the  streamwise  and  wall 
normal  Reynolds  stress  components.  Comparing 
figures  10  and  11,  it  becomes  clear  that  the  tur¬ 
bulence  model  used  in  the  computation  is  indeed 
able  to  predict  the  anisotropy  in  the  turbulent  wall 
jet. 

To  summarize,  very  good  agreement  was  achieved 
between  experimental  findings  and  the  computa¬ 
tional  results.  However,  some  of  the  shortcomings 
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Figure  6:  Normalized  skin  friction  coefficient  as  a 

function  of  normalized  downstream  distance.  - 

curve  fit  by  [18].  •  •  •  computational  results 


i r  '  ‘  ’  10* 

Re. 

Figure  7:  Local  skin  friction  coefficient  Cf  = 
rwH\/2plt)  as  a  function  of  local  Reynolds  num¬ 
ber  R(  ru-l'mVm/v. - computational  results,  •  •  • 

curve'  fit  by  [19].  -  -  -  curve  fit  by  [20] 


of  the  turbulence  model  become  apparent.  First,  al¬ 
though  anisotropic  effects  are  included  in  the  model, 
the  distribution  of  the  streamwise  Reynolds  stress 
/?n  shows  considerable  deviations  in  the  near  wall 
region.  And  second,  the  non-equilibrium  nature  of 
the  turbulent  wall  jet  cannot  be  captured  with  an 
ASM  model  of  the  form  of  equation  (3)  (or  any  other 
equilibrium  k  —  e  model). 


Figure  8:  Reynolds  shear  stress,  i?i9,  in  outer  coor¬ 
dinates.  - computational  results,  o  experimen¬ 

tal  results  at  x/b  =  70,  •  at  x/b  =  150  [14] 
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Figure  9:  Reynolds  shear  stress.  R 12,  in  wall  coordi¬ 
nates.  -  computational  results,  o  experimental 

results  at  x/b  =  70,  •  at  x/b  =  150  [14] 

3.2  Effect  of  forcing  on  the  turbulent 
wall  jet 

Using  the  Flow  Simulation  Methodology  (FSM)  in 
the  unsteady  RANS  limit,  the  effect  of  large  coher¬ 
ent  structures  on  the  turbulent  wall  jet  was  inves¬ 
tigated.  Even  though  FSM  was  developed  to  be 
used  as  a  subgrid  scale  turbulence  model,  it  is  of 
particular  importance  that  computations  in  the  un¬ 
steady  RANS  limit  yield  accurate  results.  The  dis¬ 
turbances  are  introduced  into  the  turbulent  base 
flow  using  the  volume  forcing  method  described  in 
section  2.3. 

In  experiments  with  periodic  forcing  by  Katz 
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Figure  10:  Streamwise  component  of  the  Reynolds 
stress.  R\  i .  in  outer  coordinates.  — —  computa¬ 
tional  results,  o  experimental  results  at  x/b  =  70, 
•  at  x/b  =  150  [14] 


0  03 


c.  1  2  3 

y/y , 

Figure  11:  Normal  component  of  the  Reynolds 
stress.  7?jj .  in  outer  coordinates.  -  computa¬ 

tional  results,  o  experimental  results  at  x/b  —  70, 
•  at  x/b  =  150  [14] 


rt  al.  [G] .  a  skin  friction  reduction  on  the  order 
of  10 {Z  was  achieved  with  a  forcing  amplitude  of 
only  5lZ.  From  the  large  coherent  structures  ob¬ 
served  in  a  flow  visualization  of  the  turbulent  wall 
jet.  they  concluded  that  an  instability  of  the  turbu¬ 
lent  wall  jet  exists.  They  also  showed  that  the  natu¬ 
ral  unforced  turbulent  wall  jet  has  a  strong  spanwise 
coherence,  which  is  increased  considerably  by  peri¬ 
odic  forcing.  If  these  coherent  structures  are  in  fact 
nearly  two-dimensional,  a  two-dimensional  simula¬ 
tion  should  be  sufficient  to  investigate  the  effect  of 
periodic  forcing  on  the  large  coherent  structures, 


Figure  12:  Effect  of  large  scale  structures  on 

the  mean  u-velocity  profile.  -  Base  Flow. 

■  •  •  Case02, - Case03 

Case  f/Uj 
Cased  5  x  10“3 
Case02  1  x  10~2 
Case03  2  x  10"2 


Table  1:  Disturbance  amplitudes.  Forcing  fre¬ 
quency  0j  =  0.0028. 


and  as  a  consequency,  on  the  mean  flow  of  the  tur¬ 
bulent  wall  jet. 

The  forcing  frequency  was  chosen  at  0j  =  0.0028 
(/  =  17 Hz).  Table  1  shows  the  forcing  ampli¬ 
tudes  for  the  cases  described  in  the  folk  wing.  Fig¬ 
ure  12  compares  the  mean  flow  profile*  of  Case02 
and  Case03  with  the  profile  of  the  undisturbed  base 
flow  at  x/b  =  100,150,200.  In  Case03,  the  large 
forcing  amplitudes  cause  a  significant  mean  flow 
distortion.  Note  that  this  mean  flow  distortion  in¬ 
creases  in  streamwise  direction,  indicating  that  the 
disturbances  in  the  flow  actually  enable  an  energy 
transfer  from  the  mean  flow  to  the  large,  coherent 
structures. 

The  reduction  of  the  skin  friction  is  shown  in  fig¬ 
ure  13,  where  Cf  —  Tw/{\/2pU^n)  is  plotted  as  a 
function  of  downstream  distance  x/b.  As  the  forc¬ 
ing  amplitude  is  increased,  the  deviation  from  the 
undisturbed  base  flow  increases  significantly,  not 
only  in  magnitude  but  also  in  streamwise  extent. 
For  Case03,  the  skin  friction  reduction  is  compara¬ 
ble  to  the  experimentally  observed  values. 

To  shed  some  light  on  how  the  large,  coherent 
structures  influence  the  mean  flow,  figure  14  shows 
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Figure  13:  Skin  friction  Cj  as  a  function  of  down¬ 
stream  distance.  - Base  flow,  •  •  •  CaseOl, - 

Case02. - CaseB03 


the  Reynolds  shear  stress  contribution  of  the  coher¬ 
ent  motion.  uVC!  the  Reynolds  stress  of  the  mod¬ 
elled  small  scales,  i?i2,  and  the  sum  u'v’c  +  Rn 
for  Case03.  The  figure  clearly  shows  how  the  large 
scale*  structures  increase  the  Reynolds  stress  from 
very  dose*  to  the  wall  to  about  the  wall  jet  half 
width  and  have  very  little  effect  in  the  outer  part 
of  t he  shear  layer.  As  can  be  seen  from  the  figure, 
this  also  means  that  the  y  location  of  the  zero  of 
the  Reynolds  stress  is  now  displaced  from  the  lo¬ 
cation  of  the  velocity  maximum  towards  the  wall. 
The  displacement  is  in  surprisingly  good  agreement 
with  experimental  findings  (t//e12= o  ~  O-byi/2  at 
x/b  =  200).  This  shift  has  so  far  been  attributed 
to  noneqiiilibrium  effects  and  its  origin  and  impor¬ 
tance'  has  been  the  subject  of  considerable  contro¬ 
versy.  The  simulations  presented  here  indicate  that 
the  large*  unsteady  shear  layer  structures,  in  con¬ 
junction  with  their  nonlinearly  created  harmonics, 
arc*  a  major  contributor  to  this  effect.  Figure  15 
shows  a  comparison  of  R12  from  the  ASM  model 
without  large  scale  structures,  the  sum  uVc  *f  R 12, 
and  experimental  results  of  Eriksson  et  al  [14]. 

The  local  increase  in  the  Reynolds  shear  stress 
leads  to  an  increase  in  the  spreading  of  the  turbulent 
wall  jet  and  a  reduction  of  the  maximum  local  jet 
velocity  (see  figure  12).  Both  these  effects  lead  to 
the  reduction  of  the  slope  of  the  velocity  profile  at 
the  wall,  reducing  the  skin  friction. 


10s 


Figure  14:  Reynolds  stresss  R\o.  - random  part 

(modelled),  •  •  •  coherent  part, - sum  of  coherent 

and  random  parts 


Figure  15:  Reynolds  stresss  R\o  in  inner  coordi¬ 
nates.  -  Mean  flow  without  large  structures. 

•  •  ■  Mean  flow  with  large  structures,  o  experimental 
results  at  x/b  =  70,  •  at  x/b  —  150  [14] 


4  Concluding  Remarks 

In  this  paper,  the  effect  of  large,  coherent  structures 
on  the  mean  flow  of  the  turbulent  wall  jet  is  investi¬ 
gated  with  our  new  Flow  Simulation  Methodology 
(FSM).  For  these  investigations,  FSM  was  used  in 
the  limit  of  unsteady  RANS. 

Excellent  agreement  between  experimental  and 
computational  results  was  achieved.  The  results 
show  that  the  large  scale  motion  plays  a  key  role 
in  the  transport  of  momentum  in  the  wall  normal 
direction.  This  is  most  clearly  shown  by  the  dis¬ 
placement  of  the  location  of  zero  Reynolds  shear 
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stress  from  the  location  of  maximum  velocity.  It  is 
interesting  to  note  that  the  effect  of  the  structures 
in  the  shear  layer  is  felt  very  close  to  the  wall,  in¬ 
dicating  a  strong  coupling  between  the  boundary 
layer  and  the  shear  layer  regions.  The  computa¬ 
tional  results  presented  in  this  paper  suggest  that 
these  structures  are  not  just  being  convected  down¬ 
stream  passively,  but  instead  continuously  interact 
with  the  mean  flow,  extracting  kinetic  energy  and 
redistributing  momentum  in  the  flow.  In  addition, 
the  computations  presented  in  this  paper  show  that 
the  Flow  Simulation  Methodology  in  the  limit  of  un¬ 
steady  BANS  is  indeed  capable  of  accurately  repro¬ 
ducing  the  large  coherent  motion  in  the  turbulent 
wall  jet. 
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